Non-linear Microwave Impedance of short and long Josephson Junctions 
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The non- linear dependence on applied ac field (6 W ) or current of the microwave (ac) 

impedance Ru, + iX u of both short and long Josephson junctions is calculated under a variety 
of excitation conditions. The dependence on the junction width is studied, for both field symmetric 
(current anti-symmetric) and field anti-symmetric (current symmetric) excitation configurations. The 
/**\ , resistance shows step-like features every time a fluxon (soliton) enters the junction, with a corre- 

sponding phase slip seen in the reactance. For finite widths the interference of fluxons leads to some 
\Q ' interesting effects which are described. Many of these calculated results are observed in microwave 

CI impedance measurements on intrinsic and fabricated Josephson junctions in the high temperature 

superconductors, and new effects are suggested. When a dc field (bdc) or current (idc) is applied, 
interesting phase locking effects are observed in the ac impedance Z u . In particular an almost pe- 
Q ■ riodic dependence on the dc bias is seen similar to that observed in microwave experiments at very 

low dc field bias. These results are generic to all systems with a cos(0) potential in the overdamped 
^ { , limit and subjected to an ac drive. 

Of 

O . I. INTRODUCTION 

A Josephson junction (J J) displays several remarkable properties which are well known flUpfl. In the context of 
high temperature superconductors (HTS) the analysis of JJ has received new impetus for several reasons. In addition 
to the interest in fabricating barrier junctions made of HTS, there are new features observed in the electromagnetic 
response of the HTS materials which can be attributed to the occurrence of J J distributed in the material. Many of 
the microwave properties of ceramics are ascribed to the presence of so-called "weak links" whose response is JJ-like. 
Josephson junction phenomena have been observed and attributed to Josephson coupling along the c-axis between 
superconducting layers ||] , and even within the ab-plane j| . 

Much of the analysis of an elemental JJ has been motivated by quantities which are experimentally accessible. The 
best studied parameters are the zero-frequency or dc I-V characteristics, including the Shapiro steps seen in these dc 
properties when the junction is irradiated by microwaves. These analyses have been extended to study Josephson 
junction arrays Q which can be fabricated artificially using low T c materials, and are relevant in c-axis studies. 
Recently, the high frequency or microwave impedance of single crystals has shown a remarkable non-linear depen- 
t-H ' dence on the probing microwave field which can be described in terms of Josephson-like response to purely a6-plane 
currents Results such as these in single crystals, and in fabricated bicrystal thin film junctions ||, has resulted in 
the observation of new parameters, such as the dynamic or ac (microwave) resistance and reactance. These quantities 
were not addressed extensively previously because they were not of experimental interest until recently. It is therefore 
timely to carry out a detailed study of the high frequency properties of J J, and the results of such a study are reported 
in this paper. 

In this paper we present numerical simulations of the microwave impedance of short and long Josephson junctions 
under a variety of excitations or drive conditions. In the case of a short junction, the only excitation possible is an ac 
current drive. For a long junction it is possible to study symmetric (the excitation field direction is the same on both 
O ■ sides of the junction), and antisymmetric (the excitation field direction is opposite on both sides of the junction) field 
excitations, and we consider both. We also examine the response of short and long junctions in the presence of an 
applied dc field and, in the case of a long junction, the dependence of the impedance Z u on junction width d. These 
conditions are experimentally relevant and have not been analyzed in detail before. Some of the calculated results 
are consistent with experiments on crystals and films of HTS, while other new interesting effects are described which 
may be observable in future experiments. 
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II. NONLINEAR IMPEDANCE OF SHORT JUNCTION 



The usual resistivitely-shunted JJ (RS J) model typically applies to the short junction, in which the junction width 
is smaller than the Josephson penetration depth Aj. In this limit the flux in the junction is uniformly distributed. 



f 



A. ac current driven RSJ 



Recently several authors have calculated the non- linear microwave impedance of a short RSJ |7H10(], subjected to a 
pure ac drive current represented as I sin(u;£). In this model the Josephson junction is described by an ideal junction 
J shunted by a resistance R and capacitance C to form a parallel circuit. The total current /, flowing through the 
circuit has contributions from supercurrent Is, normal current In and displacement current Ijj. While the normal 
current is caused by the flow of normal electrons in the barrier while the displacement current is due to the time 
varying electric field across the junction. Therefore, for the total current / one can write 

I = I S+ I N+ In=I c M*) + |^ + (1) 

where, </> is the junction phase and I c is the critical current of the junction. $o is the unit of a flux quantum, and the 
junction plasma frequency, w p is given by \j2 r Kl c j <&$C . For low frequency drive currents (to << oj p ~ 100GHz) the 
above equation can be written as, 

(3^ 1 d<j)/dt + sin(» = i u sin(wi). (2) 

Here, (3 = 2eI c R/h, R is the shunt resistance and, i w = I /I c . It is convenient to do the analysis in dimensionless 
form 

dip/dr + sin((/>) = i u sin^r). (3) 



where r = (3t, and Q — lo/(3. The nonlinearity arises from the sin(<^>) term in Eq. 1. Eq. 1 is solved numerically for 
(j>. The ac resistance, R u and reactance X u are calculated from the relation Q 

Z u =R u +jX u =[ 7 —) <f> e> ut dt. (4) 



Fig. 1, shows the normalized ac resistance and reactance versus microwave drive current i u for different f2 = 
0.01, 0.1, 1. As can be seen from the Fig. 1(a), when Vt is small, the ac resistance shows several characteristic features: 
a threshold field and a subsequent parabolic rise, above the threshold, with step-like structures. The threshold current 
where a sudden increase in R u starts is defined as critical current i uc . The step height decreases with increasing i u . 
The reactance directly shows dynamic phase slips occurring as the drive current is increased(Fig. 1(b)). The sharp 
steps occur at values of i u for which there are bifurcations in the solution Eq. 3, i.e. 2n phase slips in the gauge- 
invariant phase difference </> across the junction. These bifurcations can cause sudden changes in V(t) during one cycle 
of the alternating current and hence result in steps |llj]. When becomes large, both the threshold and steps in R w 
disappear, and the amplitude of the phase slip in X u becomes larger. Also, the phase slip moves to larger values of 
i u and the spacing between bifurcations becomes larger which results in R u and X u approaching R (Fig. 1). 



III. NONLINEAR IMPEDANCE OF LONG JUNCTION 



A Josephson junction is considered long when the junction width is larger than the Josephson penetration depth, Aj. 
In this case, for applied magnetic field more than the junction critical field, H c u = (3> /47rAj) In (\j/ti) p2| , ^3[ , the 
field enters the junction in the form of Josephson fluxons. A,/ = y/ ^q/[2tt^l J (2Xl + tjj] is the Josephson penetration 
depth, and Al is the London penetration depth. The RSJ model, therefore, cannot give an adequate description of 
the ac or dc field responses. Fig. 2 gives a schematic representation of a long junction for both symmetric excitation, 
when the microwave field b u is applied symmetrically on either side of the junction, in which case the induced current 
is anti-symmetric, and antisymmetric excitation due to the ac field b u , whence the microwave drive current flows 
symmetrically through the junction. 

The relation between the phase difference across the junction and the magnetic field along the junction is given by 
B z (x, t) = — 27r (2A°+t .) d ^dx ,t ^ j the relation between (j) and electric field across the junction is E y (x, t) = ■^-■§[4' { x i t) 
and the current density is J y (x,t) — J c sin0 + E y /pi. Here, the first term represents supercurrent density and the 
second term normal leakage current density, J y (x, t) is related to B z and E y by Ampere's law with a displacement 
current. By combining the above relations, for frequencies smaller than the plasma frequency, the model that represents 
the junction can be written as Q. 
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[3 = 2'KpitiJ c / ' (j) is a damping parameter. Here, U is the junction thickness, pi and e are resistivity and permittivity of 
the junction material, J c is the critical current density of the junction, and 4>o is the flux quantum. In dimensionless 
form the equation is written as 

OX*' fr=^- (6) 



where X — x/Xj, and r = 

The detailed study of the long junction depends on the parameter values and on the boundary conditions, which 
will be introduced in the following for the particuler problem considered. When the the applied microwave field 
B u = B $mujt is along z-axis and junction width infinite, the boundary conditions at x — and x — d are 



dx 



-bcj sin(f2r), 



x=0 



— 

dx 



0. 



(7) 
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Where, b u — 2itB (2Xl +ti)/<j> a , fl = to/ '(3. Eq. 7 is solved numerically for (j> and, as in the case of small junction. 
Using Eq. 4 and boundary conditions the expressions for and are obtained. 

Fig. 3 shows Ruj/Rn and X u /R n versus for different values of ft = 0.01, 0.1 and 1, for an infinite long junction 
with R n = p piUj(2\L + ti)/(2ti). When f2 is small, the variation in i? w is similar to the RSJ model calculation, in 
that it has a threshold and a parabolic rise with step-like structures. However, no clear phase slips are observed in 
X u versus b^. When fl is large, at the threshold the jump becomes higher and for further increase in b u the steps 
become smooth. When is very large the steps are washed out. The ac reactance X^ also shows a similar transitions 
from a step-like response to a smooth curve. Indeed the overall behavior of i? w and X^ are similar, this is essentially 
due to the "bulk" nature of the configuration. In the long junction the ac impedance is associated with the vortex 
- antivortex propagations. For short junctions the notion of a vortex - antivortex pair does not appear because the 
Josephson vortex has a finite size of the order 2J c Xj. Therefore, the physical mechanism leading to step structures in 
short and long junctions are different. For long junctions the steps are due to synchronous coupling of the oscillating 
vortex or antivortex with the ac current fl4| . Essentially, the steps in Fig. 3 occur due to the changes in the number 
of vortices in the junction and the step height is dependent on the microwave power. For higher values of Q vortex 
entry takes place at higher field and hence steps become broader and are moved to larger values of b^. 



In the limit the junction width d ~ A,/, the boundary condition 
ways, the symmetric and antisymmetric cases (see Fig. 2). 



7^ 0. This problem can be studied in two 



A. Symmetric field excitation (antisymmetric current drive) 



When a microwave field is applied on either side of the junction ( Fig. 2a), then this corresponds to a symmetric 
microwave excitation. The boundary condition in this case is 

= -6 w sin(nr). (8) 

x=0,d 

In Fig. 4 R u /R n versus b u responses for three different values of fl and four different values of the width, d/Xj = 5, 
10, 50 and oo, are shown. When O is small ( Fig. 4a), R u is small for low values of d/Xj compared to the high values 
of the latter. R^ approaches R at high values of b u for infinite junction width. The threshold and the steps can be 
seen clearly for all values of d/Xj. On the contrary, when Vt is large (Fig. 4c), it can be seen that there is almost no 
effect on R u with the change in the junction width, except when junction width is too small. 

Both at high field (b^ — > oo) and large damping (f2 ~ l)the junction behaves like a normal metal with R u 
approaching its asymptotic value R n . The result can be understood by considering the fact that at these two 
conditions, the R^ response is dominated by loss arising from normal current than the supercurrent. The manifestation 
of dominance of normal contribution can be seen in Fig. 4(c) in two ways. Firstly, the R u response at larger damping 
is smooth sans any threshold and step like features. Secondly, since the normal current contribution to R^ depends 
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on the Josephson skin depth, which is invariant for any fixed Cl, variation in the junction width does not have much 
affect on R u , and therefore, results in the same change in the magnitude of R u at various d/Xj values. 

However, for all the values of Cl when d/Xj is too small (d/Xj = 2), it can be seen that is small. In fact, as d/Xj 
approaches zero, the junction effects should not be observed as the total current gets canceled out due to asymmetry 
on either sides of the junction. 

B. Antisymmetric field excitation (symmetric current drive) 

If an ac current is passed through the junction ( Fig. 2b), then the direction of the induced magnetic field is 
opposite on both sides of the junction making this configuration an antisymmetric field-driven case. The corresponding 
boundary conditions are 

= b LU su\(ClT). (9) 

x—d 

In Fig. 5 plots of R^/Rn versus for d/Xj — 5, 10, 20 and oo, at three different values of Cl, are shown. When 
Cl is small, compared to the infinite width case (Fig. 3), R u is larger for small widths. It is worth noting here that 
this response is quite opposite to the symmetric field excitation, discussed in the earlier section, where i? w is larger 
when the width is large. When Cl is large enough (Cl = 1), as in the symmetric case, the change in width does not 
show much affect on R u , unless the width is too small. The interesting case, here, is when Cl is in the middle region 
(CI = 0.1, Fig. 5b). In this case for junction width several times the Josephson penetration depth (d/Xj — 5), the R^ 
response shows a threshold and rises parabolically with the step-like structures for increasing b u . It is similar to the 
result of short junction RSJ calculation and long junction small Cl case. As the junction width increases a few times 
more (d/Xj = 10), R u shows sudden jumps as b u increases. This jumpiness disappears for higher values of junction 
width, and eventually, when d/Xj = oo the steps become broader and smooth. 

These results can be explained by invoking the fact that in the antisymmetric field excitation since the currents on 
both sides of the junction flow in the same direction, as d/Xj decreases current density increases, for any particular 
field. Such a current crowding causes i? w to increase for smaller values of d/Xj. Also, the difference between the step 
patterns in the symmetric and antisymmetric field excitations is due to different vortex or antivortex interactions with 
ac current, as they move across the junction width. 

IV. NONLINEAR IMPEDANCE OF DC BIASED (CURRENT OR FIELD) JUNCTION 

An important class of experiments is where the microwave absorption is measured in the presence of a dc magnetic 
field which is varied. For this reason, we consider below the case of a JJ which is biased by a dc current or magnetic 
field. Naturally, the short junction is biased by a dc current, and the infinitely long junction by a dc magnetic field ( 
we consider the case of symmetric excitation only), and we calculate the ac impedance as a function of the dc bias. 

A. Short junction : dc current biased, ac current driven 

If a dc current, Id c is superposed on i u , the driving current becomes I sm(u>t) + Idc- The corresponding junction 
equation can then be written as 

dip/dr + sm(<f>) — i u sin(Or) + idc- (10) 

where id c = Idc/Ic- In Fig- 6, we plot R u and A w versus i^ for various values of idc — 0, 0.5, 2, and 10 at fixed 
Cl = 0.1. There are two regimes of behavior noticeable, corresponding to idc > iu>c or idc < 

When idc = 0, the critical current i^, c (idc = 0) is found to be l.l(Fig. 6). When the dc drive current idc is smaller 
than the zero dc critical current (i.e. < idc < 1-1), compared to the response of R u versus i u for id c — 0, it can be 
seen that the critical current i wc decreases and that both the step size and step rising pattern are changed. In Fig. 7, 
we plot the critical current i uc as a function of idc for Cl — 0.1. A linear fit i uc = aide + b yields a w — 1 and b w 1.1. 
Note however that this critical current value is weakly frequency-dependent. 

When idc is larger than the zero dc critical current (i.e. idc > 1-1), initially apparent jumps appear and with further 
increase of the ac drive current i u , R u begins to decrease with sudden intermittent steps and eventually approaches 



— < 

dx 



-b u sin(fh-), 



—c 

dx 
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the R u {i u — > oo, idc — 0). This regime idc, iu> > i^c is very interesting as it appears to show strong phase-locking 
effects in the ac impedance, similar to the Shapiro steps seen in dc I — V characteristics. 

For very large values of idc {idc = 10), as expected, the effect of the microwave current is negligible and the ac 
resistance approaches the shunt resistance R. 

In the case of ac reactance X u , which is plotted in Fig. 6 bottom panel, the phase slip is not affected by the increase 
in the dc current idc, however, the amplitude of phase slip reduces dramatically for larger values of idc- 

In Fig. 8, we plot i? w and X u versus idc for different microwave drive currents i u = 0.5, 1, 2, 10. It can be observed 
from the figure that for all values of as idc increases, R u initially increases, reaches a maximum and approaches 
R eventually. On the other hand, X w decreases initially with the increase in id c - The occurrence of phase slip can be 
seen clearly when idc is small. Large values of either i u or id c diminishes the effect of phase slip. 

B. Long Junction : dc magnetic field biased, ac field driven 

The presence of dc magnetic field, Bd c is studied for infinite long junction because in this case the R u response is 
same in both symmetric and antisymmetric excitations, as can be seen in Figs. 4 and 5. The calculations are, however, 
done for symmetric excitaion. When Bd c is applied in addition to the microwave field, the boundary condition becomes 

= -b u sm(ut) - b dc - (11) 

where, bd c = ^Bdd^^L + U)/4'o- Considering Q = 0.1 and d/Xj = oo case, the R^/Rn versus 6 W response at various 
bdc = 0, 0.1, 1 are shown in Fig. 9. In Fig. 10 the plots of i? w /i? n versus bdc for two different b^ values, b^ = 0.1 and 
1, are shown. It can be seen in Fig. 9 that when the dc field is high, R^ at low ac field is high and as 6 W increses it 
decreases and reaches the asymptotic value R u when b u — > bdc- This high ac resistance is due to the loss arising from 
vortex flow whose response is dependent on the intrinsic pinning in the junction. The same mechanism results in i? w 
increasing beyond R n as the dc field is increased. 

A very intriguing result is the case for intermediate b u = 0.1, shown in Fig. 10, for the parameter Q, = 0.1. The 
field dependence is strikingly regular, almost periodic. As noted below, this is similar to experimental data seen in 
microwave absorption experiments on YBCO single crystals. 

V. REMARKS 

The recent observation that bulk single crystals of high temperature superconductors behave like Josephson junc- 
tions in their microwave response, and are even of higher quality than fabricated junction was the original 
motivation for this work. In Fig. 11 the experimental microwave surface resistance variation as a function of mi- 
crowave field (H r f) recorded on a high quality single crystal of Y Ba2Cu^Oj-s, at temperature of 3K , is shown. It is 
important to note that the R s variation mimics a single long junction response for symmetric excitation (see Fig. 3), 
with Q = 0.03. Since the measuring frequency is 10GHz, this gives a damping frequency of 330G-ffz equivalent to an 
I C R = 0.8 meV. A detailed comparison of experimental data with simulations will be presented elsewhere. 

Very interesting results which can be attributed to JJ's have been observed in BSCCO pB] , ^6| . There the microwave 
surface resistance shows sudden jumps or transitions between apparently quantized values when the microwave field 
H^W c-axis. However, it is probably necessary here to also include the inertial term even in the RSJ model, which 
leads to a very rich dynamics and apparently chaotic behavior. 

Another class of experiments which have been analyzed using an RSJ model is the Modulated Microwave Absorption 
experiments |l7) , in which the microwave absorption is measured in the presence of a dc bias field, along with a very 
low frequency (< 100Hz) modulation field. Here again the properties calculated were in terms of dc quantities, and 
not the true microwave impedance as is done here, and hence the analysis should be redone along the lines of the 
present paper, as we have in the previous section. 

Indeed, as Fig. 10 shows, a striking nearly periodic behavior on the microwave bias field can be observed in the 
caculations. Similar periodic behavior has often been observed in modulation experiments on single crystals ]lq| . Our 
calculations shows that this type of behavior can occur in long junctions for a judicious choice of the parameters 
and b u . 

A surprising theme that emerges from the experiments mentioned above, when compared with calculations as 
presented in this paper, is that the experimental data can often be well described by the response of a single junction, 
rather than the multitude of "weak links" one would expect due to inhomogeneities. Thus the experiments and 
analysis taken together strongly suggest the presence of a single periodic (sin(</>)) response of the superconductor to 
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applied magnetic fields. One can imagine that this arises from a single dominant weak link in the system, or of many 
weak links acting coherently. As pointed out in ref. the degree of coherence is striking and strongly suggests a 
fundamental mechanism, such as the possibility of bulk coupled superconductors. 

Previously, several authors have modeled Joscphson junction response in granular superconductors in terms of 
strongly coupled and weakly coupled Josephson junctions. Experimentally, it is found that HTS of different families 
often respond to microwave excitations in a similar manner, including ceramics and films and different materials like 
YBCO and BSCCO. Hylton et al. [^9| have used a resistively shunted junction (RSJ) model to describe the residual 
microwave losses in thin films. In this model the grains are considered to be purely inductive and the grain boundary 
weak links the RSJ's. This model has later been expanded and used by various authors [pp[-p2[. However it is now 
clear that this early analysis needs to be replaced by the more exact analysis such as presented in this paper and by 
McDonald et. al. §]. 

While the discussion in this paper has been entirely in the context of superconducting Joscphson junctions, we are 
aware that the fundamental equations are applicable to a variety of nonlinear systems. These results are generic to 
all systems with a cos(0) potential in the overdamped limit subjected to an ac drive. Hence the results presented 
here have a general validity beyond superconductivity, wherever the ac impedance is a significant quantity which is 
experimentally accessible. A particularly interesting system is that of striped metals, where the dynamical response 
is again of the Josephson type [ p3| , and where all the analysis obtained here should apply . 

In summary, we have presented the nonlinear ac resistance and reactance analysis for both short and long junctions 
under various conditions. In the case of a short junction, the higher damping leads to less pronounced nonlincarities. 
The steps observed in the ac resistance are due to 2ir phase slips in the gauge invariant phase difference across the 
junction. The application of dc field reduces the critical current and when the dc current is more than the ac current, 
phase locking effects are seen in the ac impedance, similar to Shapiro steps in the dc I — V characteristics. 

In the case of a long junction, the step-like features in ac impedance appear due to changes in the number of fluxons. 
When the damping is large, the steps are pushed to higher values of current as the vortex entry takes place at high 
currents. In the symmetric and antisymmetric filed excitations it found that, at large damping or high fields, the 
ac resistance response is dominated by the resistive (normal ) contribution and hence leads to same magnitude and 
variation of R^, for various values of the junction width. Overall, for both short and long junctions, in high current 
and large damping limit, the ac resistance approaches the limiting dc resistance value and the nonlinearities in its 
response are less pronounced. These interesting features may be observable in experiments on crystals and films of 
HTS and other superconductors. 
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VII. FIGURE CAPTIONS 



FIG. 1. RSJ model calculation of ac (a) resistance and (b) reactance versus i^ of a short junction for three different values 
of Q, the damping parameter. The step like features seen in R^ are due to 2-7T phase slips in the phase difference across the 
junction. 



FIG. 2. Cross sectional view of long Josephson junction for (a) symmetric field excitation - antisymmetric current drive, 
and (b) antisymmetric field excitation - symmetric curret drive. 



FIG. 3. (a) Ru, and (b) versus for infinite long junction for various values of Q. The step like structures are due 
to changes in the number of vortices in the junction. At higher values Q vortex nucleation takes place at higher currents and 
hence the steps are pushed to higher current values. 



FIG. 4. R^ versus 6^ for symmetric field excitation when Q — 0.01 (a) , 0.1 (6) , and 1.0 (c), and for various junction widths. 
When the damping is high and for higher currents the resistive response dominates and hence threshold step like features 
disappear in the R^ response (bottom panel). 



FIG. 5. Rlj versus b^ for antisymmetric field excitation when (a) Q = 0.01, (b) 0.1 and (c) 1.0, and for varying junction 
widths. 



FIG. 6. Ru! (a) and X u (b) versus of short junction in the presence a dc current idc- When idc is smaller than the zero dc 
critical current, i uc (i.e. < idc < 1-1), it can be seen that i c decreases. When idc is larger than the zero dc critical current 
(i.e. idc > 1-1), Rui begins to decrease with sudden intermitent steps and eventually approaches the the limiting resistive value. 
This regime idc,i^ > i^c is very interesting as it appears to show phase- locking effects in the ac impedance, similar to the 
Shapiro steps seen in dc / — V characteristics. 



FIG. 7. Variation of i wc with idc- Note that i^ c decreases linearly with slope ~ — 1. 



FIG. 8. (a) R^ and (b) X^ versus idc for various values of i^. With the increase in idc, while R^ increases initially, X t 
decreases. The occurance of phase slip is pronounced at low ac currents. 



FIG. 9. Rui dependence on b^ for an infinite long junction when Q — 0.1 in the presence of dc field. 



7 



FIG. 10. Ru, dependence on fe dc for an infinite long junction when ft — 0.1 in the presence of ac field. 



FIG. 11. Microwave surface resistance (R s ) versus microwave magnetic field (_ff r /||c — axis) recorded at 10GHz abd 3K, 
of a very high quality single of Y Bcl^CuzOt-i,. Note that the R s variation can very well be described by a single Josphson 
junction response. 
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